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We make preliminary estimates to assess whether the recently reported flyby anomaly can 
be attributed to dark matter interactions. We consider both elastic and exothermic inelastic 
\ scattering from dark matter constituents; for isotropic dark matter velocity distributions, the 

. former decrease, while the latter increase, the final flyby velocity. The fact that the observed 

' flyby velocity anomaly shows examples with both positive and negative signs, requires the 

O 

' dominance of different dark matter scattering processes along different flyby trajectories. 



The magnitude of the observed anomalies requires dark matter densities many orders of 
magnitude greater than the galactic halo density. Such a large density could result from an 
accumulation cascade, in which the solar system-bound dark matter density is much higher 
than the galactic halo density, and the earth-bound density is much higher than the solar 
system-bound density. We discuss a number of strong constraints on the hypothesis of a 
dark matter explanation for the flyby anomaly. These require dark matter to be non-self- 
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annihilating, with the dark matter scattering cross section on nucleons much larger, and the 
^ dark matter mass much lighter, than usually assumed. 
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A. Introduction 

in 
o . 

^0 . In a recent paper, Anderson et al. [if] have reported anomalous orbital energy changes, of order 

1 part in 10^, during earth flybys of various spacecraft. Some flybys show energy decreases, and 
^ ■ others energy increases, with the signs and magnitudes related to the spacecraft initial and final 

I velocity orientation with respect to the equatorial plane. Since the DAMA/LIBRA collaboration 

has recently reported an annual modulation signal interpreted as evidence for galactic halo 
dark matter, it is natural to ask whether the flyby anomalies could be attributed to dark matter 
interactions. In this paper we give some preliminary calculations directed at this question. Needless 
to say, in proceeding along this route we are assuming that the reported flyby anomalies are not 
artifacts of the orbital fltting method used in jl]. For a detailed discussion of this, and further 
references, see jsj], which concludes that the most obvious candidates for artifactual explanations 
cannot give the large effect observed. This of course does not rule out the possibility that something 
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has been overlooked, and searching for a conventional explanation of the flyby anomaly is clearly 
a line of investigation that should be vigorously pursued.^ 



B. Elastic and Inelastic Dark Matter Scattering 

Let us consider the velocity change when a spacecraft nucleon^ of mass mi ~ IGeV and initial 
velocity ui scatters from a primary dark matter particle of mass m2 and initial velocity n2, into 
an outgoing nucleon of mass mi and velocity vi, and an outgoing secondary dark matter particle 
of mass m'2 = m2 — Am and velocity V2 ■ The inelastic case corresponds to mg 7^ m2, while 
in the elastic case, = m2 and Am = 0. (The possible relevance of inelastic scattering has 
been emphasized in a recent paper of Bernabei et al. 4]; see also the book of Khlopov Q], which 
gives arguments for unstable dark matter particles and reviews proposals js] that dark matter may 
consist of "mirror" particles.) Under the assumptions, (i) both initial particles are nonrelativistic, 
so that j-uil << c, [-1121 << c, and (ii) the center of mass scattering amplitude f{9) depends only 
on the auxiliary polar angle 9 of scattering^ , a straightforward calculation shows that the outgoing 
nucleon velocity change, averaged over scattering angles, is given by 



/.^ y m2U2-m^ui ^.ui-U2 , , 

{dvi) = ■ f ht(cos6l)— — , (1) 

mi + mg 1^1 — U2I 



with t > given by taking the square root of 

9 momf, ,^ ^ ,9 Am mn 

i = 7 : TT^. TTiui-U2?+ ' 



2^2 (mi-ui + m2'a2) 



(mi + m2)(mi + m'g) mi(mi + m2)L (mi + m2)(mi + m2) 

and with (cos 6) given by 



(2) 



d9 sine cos 9\f {6)1^ 
Jo d9sm9\f{9)\^ 



^ One possibility being discussed, and raised by a referee of this paper, is that the reported anomaly may arise from 
a mismodeling of the earth's reference frame within the barycentric system, since the earth's position relative to 
the sun is not known to a precision better than a kilometer. While it will be important to test the effects of this 
imprecision on integrations of the flyby trajectory, an argument based on energy conservation suggests that it will 
be too small. The magnitude of the change in the flyby potential energy per unit mass in the sun's gravitational 
field is ~ GMqAR/A^, with AR ~ 1.4 x lO^km the distance travelled by the flyby between ingoing and outgoing 
asymptotes, with A ~ 1.5 x lO^km the earth-sun distance, and with GMq ~ 1.3 x lO^^km^s"^. The error in this 
potential energy change arising from an uncertainty 5A ~ 1km in A is then ~ 2GMqARSA/A^ ~ 10~*km^s~^. 
However, the magnitude of the anomaly in the flyby kinetic energy per unit mass is «/ • 5vf ~ 10~^(10 kms~^)'^ ~ 
10~''km'^s~^, which is four orders of magnitude larger than the error in the sun's potential energy arising from the 
uncertainty in the earth's position. A similar calculation for the moon shows that the ~ 2cm uncertainty in its 
position relative to the earth makes a contribution eight orders of magnitude smaller than the flyby anomaly. 

^ Dark matter scattering from electrons would also be expected, with subsequent sharing of the momentum change 
with nucleons, but since this is harder to model, we ignore it for the purpose of making order of magnitude 
estimates. 

^ Here 6 is the kinematically free angle between vi — {miui -\- m2U2) / {rni +7712) and ui — {miui + 77121*2)/ [mi + 7712), 
which when 7772 7^ 77*2 is not the same as the angle between the incident and outgoing nucleon in the center of 
mass frame. 
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In the elastic scattering case, with Am = 0, = m2, these equations simpUfy to 

{6vi) = -2 {ui-U2){sm\e/2)) . (4) 

mi + m2 

In the inelastic case, assuming that Am/m2 and m'2/m2 are both of order unity, the equations are 
well approximated by 

(OTl) ~ 7-;: ^TT 7 St c(cOS^) . (5) 

pi — \ mi^nii + mgj / 

Since ui and U2 are typically of order 10 km s^^, the velocity change in the inelastic case is larger 
than that in the elastic case by a factor ~ c/juij ~ 10'^. 

To get the force per unit spacecraft mass resulting from dark matter scatters, that, is the 
acceleration, one multiplies the velocity change in a single scatter {5vi) by the number of scatters 
per unit time. This latter is given by the flux \ui — U2\, times the scattering cross section a, 
times the dark matter spatial and velocity distribution p(^x,U2)- Integrating out the dark matter 
velocity, one thus gets for the force acting at the point x{t) on the spacecraft trajectory with 
velocity ui = dx(t)/dt, 

6F = j (fu2{5vi)\ui - U2\(yp{x,U2) ■ (6) 

Equating the work per unit spacecraft mass along a trajectory from ti to tf to the change in kinetic 
energy per unit mass (assuming that the initial and final times are in the asymptotic region where 
the potential energy can be neglected) we get 

d-{vf - vf) =Vf ■Svf= r dt{dx/dt) ■ 6F 
2 Jti 

*/ /■„ 

dt d U2{dx/dt) ■ {5vi)\ui — U2\(Tp[x,U2) 

(7) 

To get the vectorial change in velocity is more difficult; one must solve the perturbed orbital 
differential equation (taking here the center of the earth as the origin of coordinates), 

d^Sx GM^f 3x^6xx\ 

— \^) ■ 

One can check that taking the inner product of this equation with dx/dt, integrating over time, and 
integrating by parts twice, again gives the energy conservation relation Vf ■ 5vf = j^* dtdx/dt ■ 6F. 

Consider now the case of a dark matter density that has an inversion invariant velocity distri- 
bution, so that p(x,U2) = p{x, — 7x2)- From ([4]) and ([7]), we see that in the elastic case, the fiux 
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weighting factor favors U2 being oppositely directed to ui = dx{t)/dt, and so the flyby velocity 
change, integrated over the dark matter velocity distribution, is oppositely directed to dx{t)/dt. 
Hence, as expected for elastic scattering, one gets a positive drag coefficient and the net effect is a 
reduction in spacecraft velocity. Turning to the inelastic case, where the flux factor in ([7|) cancels 
the denominator \ui — U2\ in (l5|), the integration over U2 leaves only the term ui = dx(t)/dt, and 
so in this case the flyby velocity change, integrated over the dark matter velocity distribution, is 
parallel to dx{t)/dt when {cos 9) > 0. So for forward dominated exothermic inelastic scattering, 
the drag coefficient is negative and the net effect is an increase in spacecraft velocity, while for 
backward dominated inelastic scattering, the drag coefficient is positive, as in the elastic case. 
Since the observations reported in l| show cases of increased velocity, and of decreased velocity, a 
dark matter explanation (assuming an approximately isotropic velocity distribution) requires the 
presence, in differing proportions on different trajectories, of inelastic forward dominated scatter- 
ing, and of either elastic or inelastic backward dominated scattering. This could be achieved in a 
two-component dark matter model, with differing spatial densities p(x, U2) governing the inelastic 
and elastic scatterers. Another possibility is a single dark matter component with an anisotropic 
velocity distribution, undergoing inelastic scattering, and possibly also elastic scattering as well. 
Detailed modelling will be needed to see which possibilities are viable.'^ 



C. Quantitative estimates 

Let us now turn to some quantitative estimates. To get a velocity change of order 10^^ of the 
spacecraft velocity over a time interval T one needs 

W-^ ^Tfap\{6vi)\/\vf\ , (9) 

with / the average flux, p the average dark matter density, a the scattering cross section, and 
K^iTi)! the magnitude of the single scattering velocity changes given, in the elastic and inelastic 
cases, by ^ and ^ respectively. This gives an estimate of the required product of mean dark 
matter density times interaction cross section, 

apr^W-^\vf\/{Tf\{5vi)\) . (10) 

* In the inelastic case, the early universe populations of the dark matter primary of mass m2 and the dark matter 
secondary of mass m'2 will be in an equilibrium resulting from dark matter scattering from nucleons or quarks. 
Hence it is not unreasonable to assume that populations of both types of particles could survive to the present 
epoch, as would be needed for a two-component model. 
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Anderson et al. 



ll report that for the NEAR spacecraft flyby, the velocity change occurs during 
an interval T = 3.7h ~ lO^s when the spacecraft could not be tracked during near earth approach. 
Taking this estimate for T and taking the mean flux as / ~ 10km s^^ = 10^ cm s^^, (jlOp becomes 

ap ^ 10-^^cm-^\vf\/\{6vi)\ . (11) 

Defining the mean dark matter mass density as pm = 1^2 P, using ^ gives for the elastic case 

ap.m ~ 10-^*^cm-i(mi + ma) > 10-^*^(GeV/c2)cm-^ , (12) 

while using ([5]) gives for the inelastic case^ 

a{Amm'2)^/^p ~ ap„ ~ 10-2°cm-^[mi(mi + m^)]^/^ > io-20(GeV/c2)cm-^ . (13) 

To estimate dark matter densities from these bounds, we must assume a value for the scattering 
cross section. For a cross section of order 1 picobarn =10"^^ cm^, we get dark matter mass densities 
Pm ~ 10^'^(GeV/c^) cm~^ in the elastic case, and pm ~ 10^^(GeV/c^) cm"'^ in the inelastic case. 
For a cross section of order 1 millibarn =10~^''cm^, which would require dark matter masses much 
below a GeV, we get corresponding dark matter mass densities pm ~ 10^^(GeV/c^) cm~^ in the 
elastic case, and pm ~ 10''(GeV/c'^) cm~^ in the inelastic case. These dark matter mass density 
bounds are orders of magnitudes larger than the estimated galactic halo dark matter mass density 
of 0.3(GeV/c^)cm~'^, but are still many orders of magnitude smaller than the earth mass density 
of about 3 X 10^^(GeV/c2)cm-^ 

Can such large dark matter densities exist in orbit around the earth? In a separate note 8||, 
we have pointed out that by comparing the total mass (in gravitational units) of the earth-moon 
system, as determined by lunar laser ranging, with the sum of the lunar mass as independently 
determined by its gravitational action on satellites or asteroids, and the earth mass as determined 
by the LAGEOS geodetic survey satellite, one can get a direct measure of the mass of earth-bound 
dark matter lying between the radius of the moon's orbit and the geodetic satellite orbit. Current 
data show that the mass of such earth-bound dark matter must be less than 4 x 10~^ of the earth's 
mass, giving an upper dark matter mass limit of 1.3 x 10^'^GeV/c^. To explain the flyby anomalies, 
earth-bound dark matter would have to be concentrated within a radius of about 70,000 km around 



For an inelastic exothermic reaction, the cross section increases as 1/u for smaU incident velocities v [?!, and thus 
the product av is what is well-defined near threshold. Rewriting (|13p in terms of af, with / ~ 10km s~^ the flux, 
we have cr/p„ > 10~"(GeV/c2) s'V 
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earth, which contains a volume of ~ 1.4 x lO^^cm^; for the dark matter mass within this volume 
not to exceed 4 x IO^^Mq, the mean dark matter density would have to be bounded by about 
10^'^(GeV/c^)cm~'^. By the above estimates, this would correspond, in the inelastic case, to a 
cross section a > lO^^^cm^, and in the elastic case, to a cross section a > 10^^^ cm^. These cross 
sections are much larger than usually assumed for the interactions of dark matter with nucleons, 
but can be compatible with existing bounds on dark matter interaction cross sections if the dark 
matter mass is much below a GeV. 



D. Accumulation cascade 

Because earth-bound dark matter mass densities of order pm ~ 10'''(GeV/c^)cm~^ or larger 
greatly exceed the estimated galactic halo dark matter mass density, a mechanism for concentrating 
dark matter near earth would be needed. One possibility is an accumulation cascade, in which 
solar system-bound dark matter is accumulated over the lifetime of the solar system, and then 
this enhanced dark matter density leads to a further accumulation near earth. Bearing in mind 
that it is an open question whether there are efficient mechanisms for dark matter capture by the 
solar system or earth Q], [2], [h] we nonetheless proceed to estimate whether such a mechanism, 
with a high capture fraction, could lead to the dark matter densities needed to explain the flyby 
discrepancies. We note also that Frere, Ling, and Vertongen [13] have pointed out that local dark 



matter concentrations in the galaxy may have played a role in the formation of the solar system, 
which could give another mechanism for producing a higher sun-bound or earth-bound dark matter 
density than the mean galactic halo density. 

Let us start with the solar system, which is moving through the galaxy at a velocity of 
Vg.s. ~ 220km s^^, with the local galactic halo dark matter approximated by a Maxwellian velocity 
distribution with a r.m.s. velocity of similar magnitude. Let /s.g. be the probability of capture of 
a dark matter particle near a solar system earth orbit of radius A = la.u. ~ 1.5 X lO^km. Then 
assuming particles captured in an annulus of radius A and area 2TTAdA over the solar system life- 
time Tg.s. 1.5 X 10^'^s are redistributed, over time, into a volume AirA^dA, the captured particle 
mass density at radius A would be 

Pm;s.s./Pm;hido ~ ^^Vg.s.Ts.s. ~ 10^^/s.s. • (14) 

So for /s.s. of unity, a very large concentration of dark matter particles in the solar system would be 
possible. In fact, the known limits on a local excess of solar system dark matter [3], [3] are about 
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3 X 10^ times the galactic halo mass density, so fs.s. in ()14p could be of order 10~^ at most. We 
remark in passing that an enhanced solar system density of dark matter particles would show up 
as a daily sidereal time modulation of dark matter particle counting rates in sufficiently sensitive 
experiments of the DAMA/LIBRA type, just as the galactic halo dark matter density is detected 
by DAMA/LIBRA as an annual modulation [14] in the counting rate. 

Given an enhanced solar system dark matter density, we can now make a similar estimate of the 
maximum possible capture density in an earth orbit, by replacing /s,s by the corresponding earth 
capture fraction /e, replacing Vg.s. by the orbital velocity of earth around the sun Ve ~ 30km s~^, 
and replacing A by the earth orbit radius relevant for the flyby anomalies, ~ 7 x lO'^km. This 
gives 

/e 

Pm;e/Pm;s.s. ~ 4^^e7's.s. ~ 2 X , (15) 

where we have divided by an extra factor of 2 since we are assuming that the solar system dark 
matter density is linearly increasing over its lifetime. So if the solar system dark matter density 
were equal to its upper bound, and fe were of order unity, the earth-bound dark matter density at 
or below the radius relevant for the flyby anomalies could be as large as ~ 10^^ times the galactic 
halo density. So even with small values of fe, one could attain large enough values of dark matter 
density to explain the flyby discrepancy if the interaction cross section were large enough. 



E. Constraints 



In addition to having to provide a large enough dark matter density, such a mechanism would 
have to lead to a dark matter spatial distribution satisfying significant constraints. We shall 
consider three types of constraints, (1) constraints coming from data on closed orbits of satellites, 
the moon, and the earth, (2) constraints coming from stellar dynamics, and (3) constraints coming 
from earth and satellite heating. 



1. Closed orbit constraints 



We begin with an analysis of closed orbit constraints, by asking what is the most general form 
of a drag force that gives zero cumulative drag for all closed satellite orbits. Let us rewrite ([7]) for 
the work per unit spacecraft mass as 

5W = J dt{dx/dt) = J d0{dx/de) ■ 5F , (16) 
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with 9 the angle in the orbital plane between the orbit semi- major axis and the vector from the 
earth's center to the satellite, and let us define the "drag function" D{x,v = dx/dt) as 

D{x,v) = {dx/de)-5F . (17) 

Then the condition for vanishing cumulative drag over the orbit becomes 

/ deD{x{e),v{9)) =Q . (18) 

JO 

Since each pair x, v is Cauchy data that corresponds to a distinct orbit, a general solution to (jlSp 
is 

oo 

D{x,v) = ^{aismW + hiCosie) , (19) 
1=1 

with determined by x^v and with the coefficients a^, bi functions of the five orbit constants 
of motion (angular momentum vector, energy, and semi-major axis orientation) that in turn can 
be computed as functions of x, v.^ That is, (jlSp is satisfied by requiring that the Fourier series 
expansion in 9 of the drag function has no constant term 6o- For a hyperbolic orbit such as the 
flyby orbits, the cumulative energy change per unit spacecraft mass is obtained by integrating (jl6p 
from —6d to 6d^ with 26d the flyby deflection angle, giving 

1 °° 6 

6-{v] -vl)=2hoeD + 2Y,j^^^^0D , (20) 

1=1 

where we have included the possibility of a nonzero 69 . Details of how the near-earth environment 
(such as a hypothetical dark matter distribution) influence the flyby energy change appear through 
the coefficient functions h^. In particular, the fitting formula given in [1] would have to arise this 
way, through the dynamics determining the coefficients 6^, and not through the kinematics of 
requiring vanishing drag anomaly for closed orbits, corresponding to vanishing hQ? 

This analysis suggests that if the flyby effect is confirmed, there likely will be analogous drag 
anomalies in high-lying satellite orbits, since a vanishing 60 would require a "fine-tuning" in the 



^ This statement applies to x, v values that do not correspond to earth-intersecting orbits. 

^ For example, one might try a dynamical model in which there are inelastic and elastic scatterers with roughly 
similar values of scattering cross section times density, with the inelastic scatterers distributed in a prolate ellipsoid, 
elongated towards the poles, and the elastic scatterers in a oblate ellipsoid, somewhat elongated towards the 
equator. For circular satellite orbits in the overlap region of the two distributions, the positive and negative drag 
effects would cancel; for a flyby deflected from a small to a large angle with respect to the equator, the negative 
drag effects would predominate, giving a velocity increase. For a flyby deflected from a moderate to a smaller angle 
with respect to the equator, the positive drag effects would dominate, giving a velocity decrease. Such a model 
would not reproduce the fitting formula given in p|, but with appropriate shapes of the density profiles might be 
able to accommodate the six flyby data sets used to generate that flt. Clearly, this is but one example of many 
possible scenarios. 
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drag law, with cancelling negative and positive drag contributions around closed orbits. Since 
normal satellite atmospheric drag effects are proportional to the cross-sectional area of the satellite, 
whereas dark matter scattering drag (of either sign) is proportional to the mass of the satellite, it 
would be helpful to have an analysis of drag effects in existing satellites, assuming the presence 
of both area-proportional and mass-proportional components. The aim would be to see if there is 
any evidence for small mass-proportional drag contributions, or at least to place bounds on such 
contributions for use as constraints on dark matter model fits to the flyby data. 

Ignoring now the possibility of fine-tuning of the drag force that could give cancellation between 
negative and positive contributions over closed orbits, let us analyze several constraints that come 
from observation of the rate at which the radius of an orbiting body increases or decreases, which 
can be used to bound a drag force acting on it as follows. For definiteness let us consider the 
earth's orbit around the sun, since other cases can be obtained from this by appropriate substi- 
tutions. Approximating the earth's orbit as circular, the total energy (potential plus kinetic) is 
E = —GM,^Mq/{2A), and the orbital velocity is given by fg = GMq/A, from which one easily 
derives that 

dE 1 dA 



(21) 



Letting dA be the change in A over a single orbit, the left hand side of ()2ip is given by 2t:A times 
the force per unit mass-energy, which by use of ([5]) for the inelastic case and ^ is given by 

2'KA{m2/mi){ve/c)aps.s. , (22) 

with ps,s. the mean density of sun-bound dark matter along the earth's orbit. Again writing 
Pm;s.s. = "i2/0s.s.) from (I2ip and ()22p we get the relation 

1 dAve , . 

0"/Om;s.s. = "j 72"— "^1 • (23) 

47r yl^ c 

Taking for dA the uncertainty in the change in A over one orbit, this gives a bound on the product 
apm acting over the orbit. 

If an inelastic earth-bound dark matter scattering mechanism is responsible for the negative drag 
flyby anomalies, then the earth's motion through sun-bound dark matter will produce acceleration 
anomalies in the earth's orbit, which (assuming no fine-tuning cancellations) can be used to place 
a bound on the density of sun-bound dark matter. Current bounds on the yearly change in the 
radius A ~ 1.5 x lO^km of the earth's orbit are around 1.5 m per century [ij], or dA ~ 1.5cm per 
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orbit. Substituting this, the value of A, and Ve ~ 30kms ^ into ([23|) gives the bound^ 

crpm;s.s. < 0.5 X 10"^^ (GeV/c2)cm~^ . (24) 

However, this formula must be corrected to take into account the fact that for cross sections 
a > 10^^'^ cm^, the earth diameter exceeds the optical depth for dark matter scattering on nucleons, 
and so not all nucleons in the earth have an equal probability of undergoing a dark matter scattering. 
Letting denote the participating fraction of earth nucleons, (I24p must be modified to read 



<jpm;s.s.Fe < 0.5 X 10"^^ (GeV/c^jcm"^ . (25) 

In terms of the density of nucleons in earth /Oearth ~ 3.3 x lO^^cm^^ and the earth diameter 

-Cearth ~ 1.3 X 10^ cm, , an estimate of F is 

1 0.2 X 10"33cni2 
Fe . (26) 

Pearth-'-^earthC 

When substituted into ([25|) . this gives the bound 

Pm;s.s. < 2 X 102(GeV/c2)cm-3 . (27) 

This bound (which we emphasize depends on the hypothesis of an inelastic dark matter collision 
explanation for the flyby anomaly, and assumes no cancellation of negative and positive drag effects 
for the earth orbit) is considerably lower than the current limit of ~ 10^(GeV/c'^)cm~'^ on excess 
solar system dark matter. Taking the ratio of (|13p . which refers to pm;e, to (j27p . and comparing 
with (jl5p , we learn that the earth capture fraction in the cascade scenario must obey the constraint 

^^^02x^0:35^ . (28) 

(7 

For a = 10~33cm^, this requires the relatively large earth capture fraction /e > 0.2 x 10~^, but for 
larger values of a the requirement on /e becomes less stringent; for example, for a = 10^^ ''cm^, 
(|28D becomes fe > 0.2 x 10^^. 

We next apply (j23p to the moon's orbit around the earth. Lunar ranging [3] has established 
the position of the moon to within a post-fit residual accuracy of about 2 cm relative to an earth- 
moon distance of Am = 384, 000km. The moon is found to be receding from the earth at a rate of 
3.8cm per year, or 0.28 cm per orbit, which is explained by the action of tidal effects. Estimating 
the uncertainty in this as dAm ~ 0.07 cm for one orbit, and substituting this, the moon's orbital 



Since for an exothermic reaction the cross section varies inversely with velocity for small velocity, the effective 
cross section a entering (|24[) through (|26[l may be a factor of 3 smaller than that relevant for the flyby anomaly. 
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velocity Vm ~ Ikms ^, and the moon's orbit radius Am into (j23p . we find that along the moon's 
orbit we must have 

apm;e < -^^^mi ~ 10-29(GeV/c2)cm-i . (29) 

However, since the radius of the moon is about 0.3 that of the earth, and the density of the moon 
is about 0.6 that of earth, for cross sections a > 6 x 10^^'^ cm^, a correction for optical depth 
~ 6-Fc ~ 10~'^'^cm^/(T is again needed. Following the reasoning of (j25p through (j27p . we end up 
with the constraint 

Pm;e < {GeY / c'^)cm-^ . (30) 

Hence the earth-bound dark matter density at the orbit of the moon would have to be many orders 
of magnitude smaller than the dark matter density within the radius of 70,000 km relevant for the 
flyby anomaly. For example, for an inelastic cross section a ~ 10~^^cm^, the dark matter density 
at the moon's orbit would have to be 10~^ of that needed to explain the flyby anomaly, while for 
an inelastic cross section of 10~^^cm^ it would have to be a factor 10~^ smaller. 

There are also low altitude constraints coming from considering satellite orbits. The satellites 
of the global positioning system have orbit radius of 26, 600 km, and geosynchronous satellites 
have orbit radii of ~ 42, 000 km, but the orbits of these satellites have not been monitored to the 
level of precision [17] of that of the LAGEOS geodetic satellite [l8|, with orbit radius of ~ 12, 300 
km. Residual accelerations of the LAGEOS satellite, believed to arise from drag effects related to 
crossings of the earth's shadow, are smaller in magnitude than 3 x lO^^^ms"^, as compared with 
the anomalous flyby accelerations ~ 10~^10^ms~^/10^s = 10~^ms^^. Thus, dark matter densities 
at the radius of the LAGEOS orbit would have to be smaller by a factor of 3 x 10^^ than at the 
orbit radii relevant for the flyby anomaly, corresponding to a constraint, in the inelastic case, 

(7pm;e < 3 X 10'^^ {GeY / c^)cm^^ . (31) 

It would clearly be of interest to have comparable anomalous acceleration limits for the higher- 
orbiting global positioning system and geosynchronous satellites, since these come closer to the 
radius 70,000 km relevant for the flyby anomaly. 

We consider finally what the comparable figure would be for Phobos, the moon which orbits 
Mars with an orbital radius of ~ 9, 400km, with an orbital period ~ 7h40m, an orbital velocity of ~ 
2.1kms~^, and an orbital radius decay of 1.8cm y~^. Application of (j29p to this case, assuming that 
the residual uncertainty in the orbital decay after taking account of tidal effects is approximately 
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1/4 of the measured value, gives an upper bound to the dark matter density at the Phobos orbit 

apm<2x 10-2^(GeV/c2)cm-^ , (32) 

which is a factor of 100 tighter than the LAGEOS bound of (1311) . Since p^a near Phobos is 
necessarily greater than the galactic halo density of 0.3(GeV/c^)cm^^, this implies that the inelastic 
cross section is bounded by o" < 7 x 10~^^cm^. Conversely, since we have inferred from the limit on 
total earth-bound dark matter that cr > 10^'^'^ cm^, we also learn that near Phobos we must have 
Pm;s.s. < 2 X 10^(GeV/c'^)cm~^, consistent with known limits on solar system-bound dark matter. 



2. Stellar (and solar) dynamics constraints 

Other possible problems raised by postulating a sun-bound dark matter density larger than 
the galactic halo density are whether the resulting dark matter accretion on the sun unacceptably 
alters our well-understood model of solar dynamics, either (i) through additional energy deposition, 
(ii) through modifications in energy transport, or (iii) by exceeding the uncertainty in the loss of 
solar mass from radiation and solar wind. 

The possible problem (i) has been discussed in detail, through a running of stellar dynamics 
codes including dark matter capture, in a recent paper of Fairbairn, Scott and Edsjo 19[|, and 
concludes that "for a spin-dependent WIMP-nucleon cross section of o" = 10^^^ cm-^, stars only 
start to change their behavior when immersed in a dark matter density of around 10^ or 10^ GeV 
cm~^", which corresponds to crpm& ~ 10^'^'^ to 10^^^(GeV/c^)cm^"^, with p^Q the dark matter 
density near the sun. To convert to a limit on Pm]s.s.^ which we have defined as the density of sun- 
bound dark matter near the earth's orbit, we should take account of the fact that the sun-bound 
dark matter density near the sun may be higher than that near the earth's orbit. Dividing by a 
factor of A/Rq = 1.5 x 10^/7 x 10^ = 214, with Rq the solar radius, as suggested by (jl4p . we can 
write the constraint coming from [19] as 

10 '^'^cm^ 

Pm s.s. < (5 to 50)(GeV/c^)cm-2 , (33) 

' cr 

which is a highly restrictive limit on the density of sun-bound dark matter. However, several 
caveats are needed. First of all, since the cross section for an exothermic reaction varies inversely 
with the velocity, and since the dark matter velocity near the sun is ~ 300km s~^, which is about 
30 times larger than the velocity ~ 10km s~^ of possible earth-bound dark matter relevant for the 
flyby anomaly, the cross section entering (j33p is 30 times smaller than the cross section relevant for 
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the flyby, making the constraint correspondingly less restrictive. Second, for cross sections larger 
than 10"'^^ cm^, the optical depth for dark matter colliding with the sun is smaller than the 
solar radius, and so dark matter with scattering cross sections from nucleons relevant for the flyby 
anomaly will not penetrate to the solar core, changing the way in which dark matter accretion 



alters solar dynamics. Finally, the limit of [19[] assumes that dark matter is self-annihilating, so 
that the total mass-energy is deposited in the sun. For non-self-annihilating dark matter, with 
the secondary dark matter particle weakly interacting so that it escapes from the sun, only the 
recoiling nucleon kinetic energy is deposited, which is at least 3 orders of magnitude smaller than 
the primary dark matter particle mass-energy for MeV or lighter dark matter (see (j37p and (jSSh 
below). For these three reasons, we expect the sun-capture constraint of ([33|) to be substantially 
weakened, and it should pose no problem ^ 

The possible problem (ii), dark matter modification of energy transport in the sun, has been dis- 



cussed by Lopes, Bertone, and Silk [20(] and Lopes, Silk, and Hansen [2]|, making use of constraints 
coming from helioseismology. They find that "in order to be effective in heat transport, WIMPs 
must have mean scattering cross section per baryon in the range of 10~^^cm^ < 0"s < lO^'^'^cm^, 
depending on the annihilation cross section and the mass of the WIMP. The transport of energy by 
WIMPs falls rapidly outside of this range." The cross sections of interest for the flyby are mainly 
larger than the upper end of this range, and correspond to a dark matter optical depth that is less 
than the solar radius, so that dark matter particles interact with nucleons before penetrating to 
the solar core. In the inelastic, exothermic scenario, dark matter particles in their first interaction 
convert to a secondary dark matter particle, that if weakly interacting escapes from the sun and 



does not contribute to energy transport. So the scenario of [20] and 211], in which dark matter 
particles contribute either to non-local or localized diffusive energy transport, is not realized, and 
does not give constraints. Additionally, recent revisions in the chemical composition of the sun 



22l | appear now to conflict with current detailed models of the inner dynamics of the sun, since 
observed sound speeds are inconsistent with the values predicted by solar models using the revised 
chemical composition. So it now appears that the helioseismology constraints may not be as tight 



as assumed in the papers [20( and 



^ I wish to thank Pat Scott for an email pointing out that the steUar dynamics constraint hkely requires that we 
assume the sun-bound dark matter to be non-self-annihilating, and for bringing the paper il9j to my attention. It 
would be interesting to know the limits analogous to those of [19|] for various cases of non-self-annihilating dark 
matter. 
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For the possible problem (hi), one can use the dark matter mass capture rate formula^'^ 

M~ap^;Q(M0/mi)(7;L/^^dm) , (34) 

with Vesc ~ 620kms~^ the escape velocity from the sun and with Vdm ~ 300kms~^ the velocity 
of dark matter near the sun. Using the bound crp^Q < 10~^^(GeV/c^)cm~^ inferred above from 
[l^ . this evaluates to M ~ 4 x lO~^^M0y~^, which is smaller than the estimated rate [15] of solar 
mass loss from radiation and solar wind, of ~ 9 x lO^^^MQy"^, and of the same order as the 
uncertainties in this rate. This estimate assumes that the entire mass-energy of the accreted dark 
matter particle is retained in the sun; in the inelastic scenario in which the dark matter secondary 
escapes, the corresponding M will be much smaller. 



3. Other astrophysical constraints 



Various astrophysical constraints on dark matter scattering cross sections from nucleons are 
reviewed in Sec. II A of Mack, Beacom and Bertone 25|]. For a dark matter mass m2 smaller than 



a GeV, these require that the dark matter scattering cross section from nucleons should be smaller 
than about 3 x 10~^^(m2C^/GeV)cm-^, which is compatible with much of the cross section range 
inferred from our analysis of the flyby anomaly. (In Sec. II B of 25|, the authors summarize direct 
detection constraints, and show that for dark matter masses below a GeV, the entire cross section 
range between 10^^ cm-^ and 10~^''cm^ is allowed.) Constraints on dark matter interactions coming 
from primordial nucleosynthesis have been discussed by Serpico and Raffelt [2ff] and by Cyburt, 
Fields, Pavlidou, and Wandelt The former paper gives constraints on the allowed leptonic 

couplings of MeV range dark matter particles. The latter studies dark matter scattering from 
baryons, when dark matter is also self-interacting, and concludes that the dark matter scattering 
cross section from nucleons must be less than about 10~^^(m2C^/GeV)cm^, which is more restrictive 



than the astrophysical constraint of 
relevant for the flyby anomaly. 



25l | but is still compatible with much of the cross section range 



4-. Earth and satellite heating constraints 

Next, we consider constraints coming from earth and satellite heating, which we shall see are 
sensitive to the value of the dark matter particle mass m2. In making these estimates, we as- 



This is eq. (11) of f^^l, which is a simplification of (2.31) of Gould [31; when the dark matter mass is much larger 
or smaller than the nucleon mass, the capture rate is smaller than this estimate. 



15 



sume that the dark matter is not self- annihilating; however, the self-annihilating case would give 
estimates similar to the first case discussed below. We shall consider bounds on the dark matter 

n 

density near the earth's surface Pm;i?® following [25] from the earth's heat flow budget. Let us 
focus on the inelastic case, and suppose that the cross section for the primary dark matter particle 
of mass 1712 to inelastically scatter on a nucleon into the secondary dark matter particle of mass m2 
is larger than 10"'^'^ cm^. In this case, the optical depth of the earth is less than one earth diameter, 
and an appreciable fraction of earth-intersecting primary dark matter particles will interact. There 
are then two limiting cases to consider. In the first case, the cross section for interaction of the 
secondary dark matter particle of mass m'2 is also large enough for this particle to be trapped 
within the earth. The kinetic energy Amc^ ~ m2C^ of this particle is then dissipated within the 
earth, and contributes to the earth's heat flow budget. The luminosity of the earth is approxi- 
mately 44TW ~ 2.8 X 1023GeVs"\ of which roughly half is accounted for by known mechanisms. 
So a.™.i„g a da* matter mass density neat eatth with a velocity of lO^cms-. a. eatth 

geometric cross section of 47r(i?0 = 6.4 x 10^ cm)-^, and including [251] a solid angle acceptance 
factor of 1/2, we get the inequality 

^/Om;iJeC^10^cms~M7r(6 X 10*^cm)2 < ^2.8 X lO^^GeVs"^ , (35) 

which gives the restrictive bound^^ 

Pm;He < 0.06(GeV/c2)cm-3 . (36) 

In the second case, the cross section for interaction of the secondary dark matter particle is very 
small, so that it escapes from the earth without interacting. In this case only the much smaller 
kinetic energy 

<5Ti ~mi(<5?Ti)V2 (37) 

of the recoiling nucleon is deposited in the earth. From ([5]), the ratio of this energy to Amc^ is of 
order 

(38) 



I wish to thank Susan Gardner and John Beacom for bringing this issue, and reference [25|, to my attention. 
Since the moon's surface heat flow is less than half that of the Earth, similar reasoning gives a bound on the dark 



matter density near the surface of the moon that is half that of (|36p . subject to the caveats that follow. 
" The DAMA/LIBRA experiment 0| sees an annual modulation signal attributed to 2 to 4 keV nucleon recoils. If 
we assume Am ^ 7712, then (|38p gives this nucleon recoil energy for a dark matter mass 1712 of 2 to 3 MeV. 
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which gives an effect that depends on the magnitude of m2. For example, for m2 of order 10 keV, 
(f38l) is of order 0.5 x 10~^, and the bound of (f36]) is altered now to 

Pm;R^ < 10^(GeV/c2)cm-3 . (39) 

Although the average effect of dark matter collisions is to alter the forward velocity of the flyby, 
there is also a random change in velocity that is averaged out in ([7]), which would show up as an 
increase in spacecraft temperature, as well as in possible localized structural disruption. Applying 
()37p in the flyby rest frame gives an estimate of the thermal energy gain by a nucleon per collision. 
Dividing this by the velocity gain by a nucleon per collision from ([I])-®, and multiplying by the 
total velocity change in the flyby (of order ^ 10~^|ni| ~ Icms^^), gives the thermal energy gain 
by a nucleon in the course of the flyby. 



Temperature gain ~ 10 ~ -mil^iJillO |ui| 

\ovi\ 2 

In the inelastic case, ([5]) gives \5vi \ ~ m2c/mi, so that ([^0]) becomes 



(40) 



Temperature gain ~ -10 ^m2\ui\c ^ 0.2°K i ^ ] . (41) 

Similarly, in the elastic case, (HD gives \6vi \ ~ m2\ui — U2\/rni, so that ([1U|) gives 

Temperature gain ~ -10~^m2|'Ui||tti — U2I ~ 10~^°K . (42) 

These imply that the dark matter mass m2 cannot be too large, or the temperature gain by the 
flyby would be noticeable; for example, from the inelastic case ()4ip we learn that the dark matter 
mass is constrained to be substantially less than a GeV. These results also suggest that sensitive 
calorimetry in high orbiting spacecraft, and perhaps even sensitive acoustic phonon detection, could 
be used to test the hypothesis that the flyby anomalies arise from earth-bound dark matter. 

If the dark matter particles are too heavy, collisions with the spacecraft nucleons will cause 
recoils energetic enough to produce structural disruption. If we require that each individual nucleon 
recoil should not produce structural changes, then we get a condition of the form 

5Ti < -^binding ! (43) 

with -Ebinding & characteristic atomic binding energy. In the inelastic case, where (p8]) with Am ~ 
1712 gives 5Ti ~ m2C^ /mi, we then get the condition 

m2c2 < (mic^Ebinding)^^^ ~ lOOkeV , (44) 
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where for sake of illustration we have taken -^binding as 10 eV. Again, we see that dark matter 
particles, if responsible for the flyby anomalies, cannot be too massive. 

In a steady-state situation, the dark matter particle capture at the earth's surface radius 6,400 
km would have to be balanced by dark matter accumulation from solar system-bound dark matter, 
at or above the radius 70,000 km relevant for the flyby anomaly. Ignoring evaporation, which should 
be taken into account in a more careful estimate, and assuming similar dark matter velocities at 
radii 6,400 km and 70,000 km, this gives as the balance condition 

Pm;s.s. f e70^ ^ p„r.R^6A^ , (45) 

which with Pm;s.s. < 2 X 10^(GeV/c^) cm~^, and using /e < 1, gives the bound 

Pm;R^ < 2.4 X 10^(GeV/c2) cm-3 . (46) 

This bound, from the steady state condition, is compatible with that of (139p obtained from the 
earth heat flow budget. We conclude that not only must the dark matter density be much smaller 
near the moon's orbit than at the radius relevant for the flyby anomaly, but it also must be 
substantially depleted near the earth's surface, to be consistent with estimates based on earth 
capture. Whether this depletion would extend to radii as large as the 30-40 thousand kilometer 
range relevant for high orbit satellites is not clear. 

5. Selection rule 

In the inelastic scattering scenario, we are postulating a dark matter primary that scatters from 
a nucleon into a lighter dark matter secondary. If the dark matter primary only interacted with 
nucleons (or quarks), one could close the nucleon line into a virtual loop, insert electromagnetic 
vertices, and deduce a rapid decay of the dark matter primary into the dark matter secondary 
plus photons. To forbid this, the dark matter primary must have interactions with other particles, 
which cancel the virtual nucleon contribution to such photon decay loops. This requires an appro- 
priate selection rule in the underlying theory of dark and ordinary matter. For example, consider 
a symmetry operation that maps the primary dark matter field into itself, the secondary dark 
matter field into minus itself, the electromagnetic field into itself, and interchanges the nucleon (or 
quark) field with a new baryonic field. Then the symmetry would forbid the decay of the dark 
matter primary into the dark matter secondary plus photons (the loop involving nucleons would be 
cancelled by the loop involving the new baryonic field) , but scattering of the dark matter primary 
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from nucleons into the dark matter secondary would be allowed, since the symmetry would only 
relate this cross section to the corresponding cross section for scattering of the dark matter primary 
from the new baryons. 



F. Summary 

To summarize, we have made a preliminary survey of whether dark matter interactions can 
explain the flyby anomaly. Our estimates do not rule out this possibility (for example, we do not 
find a requirement that /c >> 1), but the constraints are severe. To explain the cases of negative 
drag flybys, exothermic inelastic scattering of dark matter on ordinary matter is required. The 
cases of positive drag require either elastic dark matter scattering, or an asymmetric dark matter 
velocity distribution in the inelastic case. In addition, the dark matter must be confined well 
within the moon's orbit and depleted near the earth's surface, a cascade accumulation mechanism 
is required to reach the needed dark matter density, the dark matter mass must be well below a GeV, 
the dark matter interaction cross section with nucleons must be relatively high (with an inelastic 
cross section lying between around 10^'^^ cm^ and 10^^'^cm^), and dark matter must be non-self- 
annihilating and stable in the absence of free nucleons. These constraints can be compatible, since 
for dark matter masses much below a GeV, there is little information on nucleon scattering cross 
sections, Further detailed modelling will be needed to see whether the various constraints can be 
fulfilled. 15 

One could of course take the severity of the constraints as an indication that the flyby anomaly 
must be artifactual, and this may ultimately turn out to be the case. But if the anomaly is 
confirmed, and if the DAMA/LIBRA hints of light dark matter are also confirmed, then new 
physics^^ will be required, and the scenarios sketched here represent a possibility that merits 
further exploration. 



See figures 1.7-1.9 of which summarize bounds from experiments searching for muhi-GeV dark matter, and 
[3] for Hght dark matter fits to the DAMA/LIBRA signaL For dark matter self-interactions, an analysis [i^ of the 
"bullet cluster" constrains the ratio of dark matter self-interaction cross section to mass to be a jmi < 0.7cm^g~^, 
which corresponds to o" < 10~^^cm^(m2 c^/lOkeV). For a related discussion, see S. L. Adler, arXiv:0808. 28231 
Phys. Lett. B (in press). 

We also remark that for dark matter masses close to the electron mass, scattering from electrons can lead to much 
higher capture rates than scattering from nucleons with the same cross section, since the capture rate formula given 
in eq. (2.31) of Gould [i^ scales as the inverse of the mass m of the particle from which the dark matter scatters, 
at the resonance peak where the dark matter mass is equal to m. However, the dominant capture mechanism may 
be gravitational involving three-body interactions, as simulated in [ll| . 

Another possibility, corresponding to "new physics" , is that the anomaly is real and indicates that there is some- 
thing wrong with our understanding of electromagnetism [2^ or gravitation [3Qj . 
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